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THEIR CRITICAL VALUES 
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Abstract. We present three examples to illustrate that in the continuation 
of a family of normally hyperbolic C 1 manifolds, the normal hyperbolicity 
may break down as the continuation parameter approaches a critical value 
even though the corresponding generalized Lyapunov-type numbers remain 
uniformly bounded below their critical values throughout the process. In the 
first example, a C 1 manifold still exists at the critical parameter value, but it 
is no longer normally hyperbolic. In the other two examples, at the critical 
parameter value the family of C 1 manifolds converges to a nonsmooth invariant 
set, for which generalized Lyapunov-type numbers are undefined. 



1. Introduction 

In normally hyperbolic invariant manifold theory, generalized Lyapunov-type 
numbers characterize the linearized dynamics along an invariant manifold and pro- 
vide an effective way to determine the normal hyperbolicity of the invariant man- 
ifold. They are introduced by Fenichel in his seminal work on persistence and 
smoothness of normally hyperbolic invariant manifolds for flows [3]. Analogous 
concepts are also developed for maps by others (see, e.g., [5]). Roughly speaking, 
an invariant manifold is normally hyperbolic and persists under small perturbations 
if the corresponding generalized Lyapunov-type numbers are less than some critical 
values. 

A particularly interesting problem is the continuation of a family of normally 
hyperbolic invariant manifolds with respect to some parameter until the breakdown 
of normal hyperbolicity at some critical parameter value. Intuitively, one might 
expect that near the breakdown of normal hyperbolicity, some of the generalized 
Lyapunov-type numbers of these manifolds would have to become arbitrarily close 
to their critical values. However, this is in fact not the case. In [2] Chicone and Liu 
describe a scenario in which a family of normally hyperbolic limit cycles converges 
to a nonsmooth homoclinic loop and thus loses normal hyperbolicity even through 
the generalized Lyapunov-type numbers for the whole family can be made uniformly 
bounded away from their critical values. In addition, in [4 Haro and de la Llave 
study numerical continuation of invariant tori for quasi-periodic perturbations of 
the standard map and the Henon map. They observe a situation where a family 
of normally hyperbolic invariant 1-tori cannot be continued further as their stable 
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and unstable bundles converge locally when the continuation parameter approaches 
a critical value. However, the corresponding generalized Lyapunov-type numbers 
remain uniformly bounded below their critical values throughout the process. 

Compared to the vast literature on the application of normally hyperbolic in- 
variant manifolds, little attention has been paid to this counterintuitive behavior of 
generalized Lyapunov-type numbers. In fact, to the best knowledge of the author, 
the two aforementioned examples are by far the only references in the literature. 
The purpose of this paper is to provide further examples to illustrate this behavior 
of generalized Lyapunov-type numbers and to show that how much the general- 
ized Lyapunov-type numbers are less than their critical values gives no information 
about how robust the normal hyperbolicity of an invariant manifold is. In the 
following sections, we first give an overview of generalized Lyapunov-type numbers 
and then discuss three examples involving invariant manifolds for systems of ODEs. 

2. An Overview of Generalized Lyapunov-Type Numbers 

Following the presentation of Fenichel [3], we give the following definitions of 
generalized Lyapunov-type numbers. Consider a vector field defined on R n : 

(2.1) ± = f(x), 



where x G M n , F is C r for some r > 1, and ' = 4r. Let (ft be the flow of (2.1). 



Suppose M = M U dM c M n is a C r , compact, and connected manifold with 
boundary and overflowing invariant under the flow (j) 1 . Let TM denote the tangent 
bundle of M . In addition, let N denote the normal bundle of M with respect to 
the standard inner product on R n . Then, we have the splitting 

TR n | M = TM N 

with the associated orthogonal projection 

n : TR n | M — > N . 

For each p G M, define the linear operators A t (p) and B t (p) constructed from the 
linearized flow Deft as follows: 

At(p) := D4>-\p)\t p m : T p M - T - t(p) M , 

Bt(p) := iW^-'O^liV^ : JV'(p) ^ N p ■ 

Then, the generalized Lyapunov-type numbers v(p) and <j(p) are defined as follows: 

v(p) := limsup||£ t (p)||* , 

{2 ' 3) ( x r log 11^)11 

a(p) := hmsup — ,, — . N ,, . 

The number v(jp) measures the exponential of the normal contraction rate under 
D(/) t (p), and the number a(p) compares the tangential and normal contraction rates 
under D(j) t {p). In addition, both numbers are constant along the trajectory (j)~ t (p) 
for all t > 0, and u(p) = v(p') and a(p) = cr(p') for any p' in the a-limit set 
of p. The proofs of these properties can be found in [3j [6]. In [3] Fenichel has 
also introduced other generalized Lyapunov-type numbers in the same spirit as 
the above construction but for a more general situation where the normal bundle 
of M further splits into stable and unstable subbundles. However, the examples 



considered in this paper involve only the two defined by (2.3) 
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If v(p) < 1 and a(p) < - for all p G M, then the overflowing invariant manifold 
M is said to be normally hyperbolic (or more precisely, r-normally hyperbolic) and 
M persists under small perturbations according to the following theorem. 

Theorem (Fenichel [3]). Let x = f(x) be a C r vector field on R n ; r > 1. Let M = 

M U dM be a C r , compact, and connected manifold with boundary and overflowing 
invariant under the flow of x = f(x). Suppose v(p) < 1 and a(p) < ^ for all 
p G M. Then for any C r vector field f e in a sufficiently small C 1 neighborhood of 
f , there is a manifold M e overflowing invariant under the flow of x = f e (x) and 
C r diffeomorphic to M. 

When M is a C r , compact, and connected boundary less manifold and invariant 
under the flow 0*, the above construction of generalized Lyapunov-type numbers 



using (2.2) and (2.3) is still applicable. Furthermore, suppose v(p) < 1 and a(p) < £ 
for all p G M. Then we can retain the notion of normal hyperbolicity and the 
accompanying persistence property even though M is boundaryless. Specifically, 



we can append an auxiliary variable £ G 1R to (2.1 ) to form an enlarged system as 
follows: 

x = f(x) , 

2.4 . , 

c = c 3 . 



Then M. := M x [—1, 1] c M n+1 is a C r , compact, and connected manifold with 
boundary and overflowing invariant under the flow of ( |2.4| ). It is straightforward 
to verify that for any (p, C) £ -M, v{p,C) = v(p) and 5(p,() = cr(p), where v{pX) 
and a(p, () are the corresponding generalized Lyapunov-type numbers computed 



using the flow of (2.4). Then Ai is r-normally hyperbolic with respect to the flow 



of (2.4), and it perturbs to a C r manifold Ai e that is overflowing invariant under 



the flow of 

x = f e (x) , 

c = c 3 , 

for any C r vector field f e that is sufficiently C 1 -close to /. Clearly, the cross- 
section of Ai e at £ = is C r diffeomorphic to M and invariant under the flow of 
x = f e (x). In the balance of this paper, we consider only boundaryless invariant 
manifolds. Furthermore, we take r = 1 and refer to a C 1 , compact, and connected 
boundaryless manifold M as being normally hyperbolic if v(p) < 1 and a(p) < 1 
for all p G M. 

3. Example 1 

Consider the following 1-parameter family of 2-dimensional systems: 

a = — f i + sin 0) a , 
(3.1) . V2 J 

= 0, 

where a G R, G M(mod27r), and the parameter /3 G R. Clearly, the circle 

T:= {(a,6>) : a = 0, G R(mod 27r)} 

is an invariant manifold for ( |3.1[ ) with any G R. In particular, T is a periodic orbit 
if f3 7^ 0, and it is a set of fixed points if /3 = 0. Let us compute the generalized 
Lyapunov-type numbers for T. 
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Let be the flow of (3.1) with the corresponding (3. Note that for any (3 G 



Ydip) — (0,# + ft) f° r an Y P — (O>0) £ T. Suppose (3 =fi 0. Then straightforward 



application of ( |2.2| ) and ( |2.3| ) yields that for any point p in the periodic orbit T, 
073 (p) = and 

^( p ) = e -^/o 2 "(H-^ = e -| ; 

which is also the Floquet multiplier of any periodic solution along T. It follows that 



T is normally hyperbolic with respect to the flow of (3.1 ) with any (3^0. However, 



for (3.1) with (3 = 0, T is not normally hyperbolic since 

MP) = e-(^ sine ) > 1 

for any p = (0, 0) G Y with G [^7r, ^7r]. Therefore, T loses normal hyperbolicity 
at /? = even though T is normally hyperbolic with the generalized Lyapunov-type 
numbers vp(p) = e~2 < 1 and <Jp{p) = < 1 for any p G T and for any (3^0. 

4. Example 2 

As mentioned earlier, in [2 Chicone and Liu describe a scenario in which a family 
of normally hyperbolic limit cycles converges to a nonsmooth homoclinic loop while 
the generalized Lyapunov-type numbers for the whole family can be made uniformly 
bounded away from the critical value 1. However, they have not provided an explicit 
formulation of the system that admits such a dynamical feature. Below we present 
a simple example based on the mechanism introduced in [2]. 

Consider a 1-parameter family of planar systems: 

x = - 1 

(4.1) 



y = x-x 3 + y(±y 2 - (\x 4 „ 

where the parameter c G [—0.1, 0]. It is clear that the origin (0, 0) is a saddle point 
for any c G [—0.1,0]. There are three critical parameter values: c\ ~ —0.060959, 
C2 — —0.060932, and C3 —0.060903, at which bifurcations occur. For —0.1 < 



c < ci, (4.1) has two unstable limit cycles Ti(c) and T2(c) as sketched in Figure 
|l(a)| At c = ci , a saddle- node bifurcation of limit cycles creates a new pair of limit 
cycles Ts(c) (stable) and T^c) (unstable) as shown in Figure [T(b)) For c\ < c < C2, 



Ts(c) shrinks in size as c increases (Figure 1(c)) and eventually converges to a 
figure-8 homoclinic loop of the saddle point (0,0) in the limit c — > . Note that 
the homoclinic loop exists only when c = C2, at which a homoclinic bifurcation 



occurs (Figure l(d)| ). For C2 < c < C3, Ts(c) bifurcates into two smaller stable limit 



cycles r 5 (c) and Tq(c) as shown in Figure l(e)| At c = C3, r 5 (c) and Tq(c) collide 



with Ti(c) and T2(c), respectively, in two simultaneous saddle- node bifurcations of 
limit cycles (Figure [l(f)[ ). Finally, only the unstable limit cycle T^c) remains for 
C3 < c < 0, as displayed in Figure [T(g)| 

Now consider the continuation of the stable limit cycle Ts(c) for c G [c\ + 5 : C2) 
with S being some small positive constant such that c\ + S < c^. Since Ts(c) is 
a limit cycle, the generalized Lyapunov-type number cr c (p) = for all p G Ts(c). 
Furthermore, the generalized Lyapunov-type number v c (p) = A c for all p G T 3 (c) 
with A c being the Floquet multiplier of any periodic solution 7(-,c) : R — > r 3 (c) 



along Ts(c). For the planar system (4.1), A c can be computed by the following 
formula (see, e.g., [1]) 

\ = ^y/o T(c) div/(7(t,c),c)dt 
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(a) -0.1 < c < ci 



(b) c = ci ~ -0.060959 




(c) Cl < C < C2 



r 4 (cL 




(d) c = c 2 ~ -0.060932 




r 4 (c) 




r 5 (c) X 



(e) C2 < C < C3 



(f) c = C3 ~ -0.060903 



c 3 < c < 



Figure 1. Sketches of the phase portrait of (4.1 ) for various values 



of c. Black dots denote fixed points, blue dotted curves are unstable 
limit cycles, and red solid curves are stable limit cycles or stable 
homoclinic orbits (in (d)). Note that the variations of the limit 
cycles are largely exaggerated. All the limit cycles are actually 
contained in a small neighborhood of the homoclinic loop. 
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where f((x,y),c) denotes the right-hand side of (4.1) and T(c) is the period of 
7(t, c). Note that 

lim A c = e d -/((°.°)> c ) = e~h-^ < e -o.4 . 

c— >c 2 

Thus, if ci + S is chosen sufficiently close to C2, the generalized Lyapunov-type 
numbers v c {p) and cr c (p) are uniformly bounded below 1 for any p G ^(c) with any 
c G [ci + 5, C2). However, the continuation of ^(c) has to cease at c = c^. 

5. Example 3 

In this example, we consider the continuation of a normally hyperbolic invariant 
torus in a 1-parameter family of 3-dimensional systems. During the continuation 
process, the torus continuously deforms. However, it contains the same a-limit 
set. Consequently, the generalized Lyapunov-type numbers of the invariant torus 
remain constant for all the parameter values for which the torus exists. 

Consider the following 1-parameter family of systems: 

a = — a + (3 2 sin 2 £ sin , 
(5.1) C=sin 2 C, 

= a , 

where a G M, ( G R(mod7r), G M(mod27r), and the parameter (3 G [0, 1]. Clearly, 



(5.1) with (3 = has a unique invariant torus 

T = {(a,C,0) :a = 0,(G R(modTr), G R(mod27r)} . 

In what follows, we will study the continuation of To for (3 > 0. 

It is evident that for any (3 G [0,1], the set of fixed points of (5.1) forms an 
invariant circle 

r = {(a,C,0) :a = 0, C = 0, 9 G R(mod 27r)} . 

Let (j)^ be the flow of ( |5.l| ) with the corresponding (3. For any (0, 0) G T , we 
have 

lim $(o,£0) = (0,0,0) er. 

Thus, to compute the generalized Lyapunov-type numbers for points in To, we only 
need to consider points in V since ^o(0> 0) = ^o(0, 0, 0) and <to(0, C, 0) = tf"o(0 5 0, 0) 
for any (O,C,0) GT . 

For any (0,0,0) G T, the linear variational equation along the solution trajectory 
^(0,0,0) = (0,0,0) is given by 

(5a\ (-100 

(5.2) IS'C = 

\so \i 00 




By integrating (5.2), we obtain 

^o(O,O,0) = 

Notice that at any (0,0,0) G T, the tangent space to To has a basis {^,^} 
and the normal space to Tq has a basis { ^ } • Then straightforward computation 
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using the definitions (2.2) and (2.3) yields z/q (0,0,0) = e _1 and (70(0,0,0) = 
for all (0, 0, 0) G T, which verifies the normal hyperbolicity of Tq with respect 
to 00. Therefore, for some f3* > 0, there exists a continuous family of C 1 tori 
{Tj3 : (3 G [0,/?*)} such that each TJg is invariant under the flow 0^ generated by 
( |5.1[ ) with the corresponding /?. 

For each /3 G [0, 1], let be the largest invariant set under the flow 0^ inside 
the domain U defined as follows: 

U := {(a,t,0) : \a\ < 2, ( e M(mod7r), G R(mod27r)} . 

In the following lemma, we state two properties of A^ that are useful in our subse- 
quent analysis of the continuation of Tp for /3 > 0. 

Lemma 5.1. For each f3 G [0, 1], A^ has the following properties: 

(1) // (a,C,6>) G A/3 and ( G [0, §], tten |a| < sin 2 C. 

(2) For each (a, £, 0) G A^ wit/i £ G [0, 7r) ; t/iere exists a #* G M(mod27r) s^c/i 
*Aa* lim 0Ua, C, 0) = (0, 0, 6>*) G I\ 

Proof. To prove property (1), we consider the time-reversed trajectories </>^ r (a, £, 0) 
for r > 0. Clearly, if £ G [0, f], then the ^-component of (/>^ r (a,C,#) is contained 
inside the interval [0, ^] for all r > 0. Let u := a — sin 2 £. Then at any (a, 0) with 
C G [0, |] and a — sin 2 (" > 0, we have that under the flow 0^" r , 

^ = - (-a + /? 2 sin 2 C sin - 2 sin 3 C cos () 
dr 

= u + sin 2 (" — (3 2 sin 2 ^ sin + 2 sin 3 ^ cos ( 
>u>0. 

It follows that u(t) — > oo at least exponentially along ^^('Uo+sin 2 ^, C #) as r ~ > 00 
if M (o) = ii > and ( G [0, §]. Next, let v := a + sin 2 C- At any (a,C,0) with 
C G [0, |] and a + sin 2 (" < 0, we have that under the flow 0^" r , 

^ = - (-a + /3 2 sin 2 C sin (9 + 2 sin 3 C cos Q 
dr 

= v — sin 2 (" — (3 2 sin 2 (" sin — 2 sin 3 ^ cos ( 
< v < 0. 

Then v(r) — > — oo at least exponentially along <j)^ 7 '(vq — sin 2 C,C,#) as r — > oo if 
v(0) = < and ( G [0, f]. Therefore, for (a, C, 0) G A^ with ( G [0, f], the 
boundedness of the a-component of ^ r (a, £, 0) for r — > oo makes it necessary that 
|a| < sin 2 C- 

Recall that ( = sin 2 ^. Consider any (aXo,6) £ with Co G (0,7r), and 
reparametrize the negative semi-orbit originated at (a, Co 5 #) by £ G (0, Co]? i- e -? 

{MC),C,^(0) : C G (0,Co]} := {<^ T (a,Co,0) : r > 0} . 

Then property (1) implies that |a/?(C)| < sm2 C f° r an y < C < f • It follows that 
the limit 

lim %(C) = ^(Co)+ / 
C^o+ J Co sin 2; 
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exists. Furthermore, if (a, £, 0) G with £ = 0, then a = and hence </>^(a, 0, 0) = 
(0, 0, 8) for all t G R. Thus, property (2) is true. □ 

Now we take an arbitrary (3 G [0,/?*) and compute the generalized Lyapunov- 



type numbers for each point in Tp C A^. In view of property (2) in Lemma 5.1 
and the obvious fact that r C Tp, this task can again be reduced to computing 
the generalized Lyapunov-type numbers for points in T only. Note that for any 



(3 > and any (0,0,0) G T, the linear variational equation of (5.1 ) along the 
solution trajectory 0^(0,0,0) = (0,0,0) is still the same as (|5.2|. Furthermore, 
property (1) implies that Tp is tangent to To along T. Thus the tangent space of 
Tp at any point in T coincides with that of To. Therefore, for any (3 G [0,/3*), we 
still have ^(0,0,0) = e~ x and ^(0,0,0) = for all (0,0,0) G T. It follows that 
vp{a, C, 0) = e _1 and erg ( a, (, 0) = for any (a, £, 0) G Tp as long as T/3 exists. 



However, the system (5.1) with (3 = 1 has no C 1 invariant torus that can be 
continued from Tp. To see this, we consider the forward-time images of the set L 
defined as follows: 

L := {(a,C,0) : \a\ < sin 2 ^ + 10, £ = G R(mod27r)} . 



Since ( = sin 2 (", by numerically integrating (5.1) with (3 = 1 starting at points on 
the boundary of T, we obtain a series of snapshots of (j)\(L) at various sections 

E c = {(a,C,0) : a G R, G R(mod27r)} 

with monotonically increasing ( converging to ir from below as t — > oo. We find 
that a fold develops on 4>\(L) in a neighborhood of (a, 0) = (0, 7r) in the limit t — > oo 
as illustrated by the snapshots of 0i(T) taken at with (" = 0.897T, 0.97T, 0.917T, 
and 0.927T shown in Figure [5] Suppose the family of C 1 tori {Tp : (3 G [0,/?*)} can 
be continued up to (3 = 1. Then the C 1 torus Ti must contain T, and Ti p| E 



must be contained in L by property (1) in Lemma 5.1 It follows that T x f] E^_ is 
a circle winding around the "cylinder" L exactly once. However, as ( — > 7r~, the 
cross-section Ti P| E^, which is contained in 0i(T) taken at E^, has to fold as (j)\(L) 
does in the limit t — > oo. This immediately contradicts the assumption that T x is 
a C 1 torus. 



Let (3 C be the infimum of the set of (3 with which (5.1) has no C 1 invariant torus 
that can be continued from Tp. The preceding analysis shows that f3 c < 1. In 
addition, it follows directly from the definition of (3 C that the family of C 1 tori 



Tp exists for all (3 G [0,/? c ). However, if (5.1) with (3 = (3 C has a C 1 invariant 
torus T/3 c , then vp c (aX,@) — e _1 and <jp c (aX,6) = for any (a,£,0) G Tp c since 
our previous computation of the generalized Lyapunov-type numbers for points in 
Tp is still applicable. Consequently, Tp c is normally hyperbolic, and there exists a 
sufficiently small S > such that the family of C 1 t ori Tp exists for all (3 G [0, (3 c +5), 



which contradicts the definition of (3 C . Thus, (5.1) with (3 = (3 C has no C 1 invariant 
torus that can be continued from Tp although Tp is normally hyperbolic and the 
corresponding generalized Lyapunov-type numbers are constant below 1 for all (3 G 

[o,/3 c ). 

Further analysis reveals that the breakdown of Tp in the limit /5 — > is due to 
the rotation of the normal bundle of the orbit 



E:= {(a,C,0) :a = 0, (e (0,tt), 6 = it} , 
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(a) Full view 




(b) Zoom view near (a, 6) = (0, tt) 

Figure 2. The images of <j>\{L) at C = 0.89tt, 0.9tt, 0.91tt, and 
0.927T, plotted in red, green, yellow, and blue, respectively. 



which is an invariant subset of for any (3 G [0, 1]. Specifically, linearizing (5.1) 
about E gives 

8a = — 8a — ft 2 sin 2 ^ 80 , 

C= sin 2 C, 
80 = 8a , 



where 8a G R and 80 G R are the variations in a and 6, respectively. Next, with 
the substitutions 8a = psina and 80 — pcosa, we obtain a system defined on the 
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(a) (3 = 0.65 (b) (3 = (3 C ~ 0.815 (c) (3 = 1 



Figure 3. The phase portraits of ( |5.3| ) with (3 = 0.65, /? = (3 C ~ 
0.815, and (3 = 1. Note that (" = and £ = 7r are identified 
and a = and a = tt are identified. The light-blue dots denote 
the fixed point (a, C) = (0,0), an d the pink dots denote the fixed 
point (a, C) = ( ^7r, 0) . In (a) and (c), the homoclinic orbits of 
(a,C) = (0,0) and (a, C) = ( £tt, 0) are plotted in blue and red, 
respectively. In (b), the two homoclinic orbits "collide" and are 
replaced by a heteroclinic orbit shown in mixed blue and red. 

torus {(a, £) : a G M(mod7r), ( G M(mod7r)} as follows: 

d = — cos a sin a — sin 2 a — (3 2 sinV cos 2 a , 

(5-3) 

C=sin 2 C, 

which describes the rotational dynamics in the normal bundle of E under the lin- 



earized flow D<j) l J{), tt). For any (3 G [0, 1], (5.3) has two fixed points (a, () = (0, 0) 



and (a, C) = ( ^tt, 0) . For (3 G [0, /? c ), each fixed point has a homoclinic orbit, along 



which the numbers of rotation in a and £ are and 1, respectively (see Figure 3(a) 
for the case with (3 = 0.65). As (3 approaches (3 C from below, the two homoclinic or- 
bits converge pointwise to each other. At (3 = (3 C c± 0.815, the two homoclinic orbits 
"collide" and are replaced by a heteroclinic orbit (see Figure [3(b)] ) . For (3 G (/? c , 1], 
the two fixed points are again connected to themselves by homoclinic orbits, but 
now the numbers of rotation in a and £ are —1 and 1, respectively, along both 



homoclinic orbits (see Figure 3(c) for the case with (3 = 1). 

The above analysis implies the following picture about the rotation of the normal 
bundle of E under the linearized flow Z}0|(O, C 71 ")- For every (3 G [0, 1], we obtain 
a subbundle Np of the normal bundle of E with its fibers oriented according to 
the values of a along the orbit whose backward-time asymptotic limit is the fixed 



point (a,C) = (0,0) in the system (5.3) with the corresponding (3. Clearly, for any 
(3 G [0, 1], the bundle Hp is invariant under D0^(O, C, tt), and its fibers satisfy 

( 5V^I(o,M = {A|:AeR}, 

which coincides with the tangent space of Y at (a, £, 0) = (0,0,7r). Furthermore, 
for any /3 G [0, 1], Afp is tangent to along the orbit E. For any (3 G [0,/? c ), the 
bundle Mp corresponds to the homoclinic orbit of the fixed point (a, Q = (0, 0) in 
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Figure 4. with (3 = 0.65. The black arrows denote unit vec- 
tors in the fibers of Afp at various locations along the orbit E, 
which is shown in green. For any f3 G [0, /3 C ), Tp = A^, and Afp is 
tangent to Tp along the orbit E. 



(5.3), and thus 

lim A/>| (0 ,c,^) = lim ^l(o,C,7r) • 

In this case, A^ is the C 1 invariant torus Tp (see Figure[4]for the case with /3 = 0.65). 
However, when /3 = /3 C ~ 0.815, the bundle Afp corresponds to the heteroclinic orbit 
shown in Figure |3(b)| in mixed blue and red and thus 

lim 7V /3c | (0 , c , 7r) = {A(-^ + ^):AG]R}. 

Consequently, Ap c loses smoothness at (a, £, 0) — (0, 0, 7r) and becomes a nonsmooth 
topological torus (see Figure |5|. Finally, for any (3 G (/3 C ,1], since the number of 
rotation in a is —1 along the corresponding homoclinic orbit of (a, C) = (0, 0), 
the bundle Afp has a — ir rotation along the orbit E, which causes A^ to fold in a 
neighborhood of (a, 6) = (0,7r) in the limit ( — > ir~ (see Figure |6] for the case with 
^ = 1). 
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Figure 5. with /3 = /3 C ~ 0.815. The black arrows denote unit 
vectors in the fibers of Np c at various locations along the orbit E, 
which is shown in green. hp c is a nonsmooth topological torus. 

[6] S. Wiggins, Normally Hyperbolic Invariant Manifolds in Dynamical Systems, vol. 105 of 
Applied Mathematical Sciences, Springer- Verlag, New York, 1994. 



Department of Mathematics, Cornell University, Ithaca, New York 14853 
E-mail address: gy260cornell.edu 



BREAKDOWN OF NORMAL HYPERBOLICITY 



13 




(b) Zoom view near (a, £, 0) = (0, 7r, 7r) 

Figure 6. with (3 = 1. The black arrows denote unit vectors 
in the fibers of A/i at various locations along the orbit E, which 
is shown in green. Notice the rotation of the bundle J\f\ shown in 
(a). A zoom view of the fold on Ai is shown in (b) with various 
cross-sections of Ai highlighted in light blue. 



